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Summary 

An expanded and improved version of a previously pub- 
lished benzene oxidation mechanism is presented and shown to 
model published experimental data fairly successfully. This 
benzene submodel is coupled to a modified version of a toluene 
oxidation submodel from the recent literature. This complete 
mechanism is shown to successfully model published experi- 
mental toluene oxidation data for a highly mixed flow reactor 
and for higher temperature ignition delay times in a shock tube. 
A comprehensive sensitivity analysis showing the most impor- 
tant reactions is presented for both the benzene and toluene 
reacting systems. The NASA Lewis toluene mechanism’s 
modeling capability is found to be equivalent to that of the 
previously published mechanism which contains a somewhat 
different benzene submodel. 


Introduction 

This report continues the investigation of aromatic hydrocar- 
bon oxidation mechanisms begun with the development of a 
benzene oxidation mechanism by Bittker (1991), The increased 
content of aromatics in today’s practical fuels makes it impor- 
tant to understand their oxidation chemistry because a large 
research effort is being devoted to the theoretical modeling of 
advanced concepts for high-speed, clean-burning aircraft en- 
gines. This work requires the development of simplified fuel 
oxidation models which realistically predict heat release rates 
and pollutant species emission concentrations. The latter task 
can only be accomplished after an understanding of the com- 
plete detailed oxidation mechanism has been obtained. Over 
many years of research, although steady progress has been 
made toward the understanding of aliphatic hydrocarbon 
oxidation (Wamatz (1984)), this has not been the case for 
aromatics. Only very recently, with the publication of the 
NASA Lewis benzene mechanism (Bittker, 1991) and a tolu- 
ene oxidation mechanism by Emdee, Brezinsky, and Glass- 
man (1992), has detailed modeling for experimental reactions 
of these fuels been reasonably successful. The latter toluene 
mechanism contains a benzene submodel which differs in 
several respects from the NASA Lewis mechanism. Both 
benzene oxidation mechanisms have been used to successfully 
compute experimental composition time profiles measured in 
a highly mixed flow reactor at approximately 1 100 K (Lovell, 
Brezinsky, and Glassman, 1988). The NASA Lewis mechanism 


also computed with reasonable success experimental ignition 
delay times measured behind a reflected shock wave for lean, 
stoichiometric, and rich benzene-oxygen-argon mixtures 
(Burcat, Snyder, and Brabbs, 1986). The complete Emdee 
mechanism computes species profiles for the toluene oxidation 
which agree with experimental data measured at 1200 K in the 
same flow reactor used for the Lovell work. 

This report first presents an improved version of the original 
NASA Lewis benzene oxidation mechanism. The capabilities 
of both mechanisms to model available experimental data are 
compared and sensitivity analysis results for the new mecha- 
nism are given. Differences between the Emdee and NASA 
Lewis mechanisms are pointed out. The NASA Lewis benzene 
model is then coupled with a slightly modified version of the 
Emdee toluene submodel. The new toluene oxidation mecha- 
nism is used to model not only the Emdee flow reactor data at 
1200 K but also the ignition delay time data for toluene- 
oxygen-argon mixtures which were also reported by Burcat, 
Snyder, and Brabbs (1986). The temperature range of the 
Burcat data is 1300 to 1600 K; the mixture concentrations 
range from a lean fuel-oxygen equivalence ratio (pof 0.331 to 
stoichiometric (<p= 1 .0). A complete sensitivity analysis is also 
performed to determine the rate-controlling reactions and to 
indicate which individual steps in the mechanism need further 
study. 

All computations were performed with the NASA Lewis 
Research Center general chemical kinetics and sensitivity 
analysis code LSENS (Radhakrishnan, 1994, Radhakrishnan 
and Bittker, 1994, and Bittker and Radhakrishnan, 1994). 
According to the law of microscopic reversibility, all chemical 
reactions are considered reversible, with the ratio of the for- 
ward to reverse rates equal to the equilibrium constant. 
Because the netrates, therefore, depend strongly on the thermo- 
dynamic data used to compute the reaction equilibrium con- 
stants, table I lists pertinent thermodynamic data over the 
temperature range of interest for several important species in 
the benzene and toluene oxidation mechanisms. The thermo- 
dynamic data base used for most species is that of the NASA 
Lewis Chemical Equilibrium Composition Code (Gordon and 
McBride, 1994 and McBride, Reno, and Gordon (1994)). 
However, the table gives data for several species not in the 
standard data base; these data were computed recently for the 
present work by Bonnie J. McBride at the NASA Lewis 
Research Center. Data for several C 5 species were also 
computed by Dr. A. Burcat at Lewis. It is also important to note 
that the thermodynamic data for benzyl alcohol and the benzyl 
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radical are the same as those used by Emdee, Brezinsky, and 
Glassman (1992) in their computations and were provided by 
Dr. K. Brezinsky (1993, Princeton University, Princeton, N.J., 
personal communication). The exact structures of these species 
are not known at this time. Because their thermodynamic 
properties had a significant effect on the computed benzene 
alcohol concentration, it was important to use the same struc- 
tural assumptions that Emdee, Brezinsky, and Glassman used 
for these molecules in their computations. 

Benzene Oxidation 

Chemical Mechanism 

The original NASA Lewis mechanism overpredicted the 
formation of phenol and cyclopentadiene. Additional con- 
sumption reactions for these species were published in the 
benzene oxidation submodel of the Emdee toluene mechanism. 
These reactions and other modifications were incorporated into 
a new NASA Lewis benzene submodel listed in table II, 
which contains the three constants in the rate equation 
kj = AjTp exp - ( Ej/RT ) for all reactions. The Emdee prod- 
ucts and rate coefficient parameters for the reaction of oxygen 
atom with the cyclopentadienyl radical (C 5 H 5 ) are now used. 
Instead of forming the intermediate product C 5 H 5 0, the reac- 
tion now forms directly the products C 4 H 5 and CO. In addition, 
two different pyrolysis reactions of the C 4 H 5 radical are used, 
as done in the Emdee mechanism. These are 

C 4 H 5 = C 2 H 3 + C 2 H 2 (50) 

and 

C 4 H 5 + M = C 4 H 4 + H + M (51) 

The rate coefficient for reaction (50) was modified from the 
Emdee expression by reducing the pre-exponential factor to 
obtain a vinyl acetylene (C 4 H 4 ) concentration profile consis- 
tent with the experimental results. The rate coefficient of reac- 
tion (51) was assumed to be the same as the expression used for 
the vinyl radical (C 2 H 3 ) pyrolysis, reaction (60) (also done in 
the Emdee mechanism). In addition, destruction reactions of 
the species C 4 H 4 and the propargyl radical (C 3 H 3 ) were added 
to the NASA mechanism. The molecular pyrolysis of C 4 H 4 to 
acetylene, suggested by Kiefer et al. (1988), was used as were 
reactions of vinyl acetylene with the phenyl radical and H-atom 
reported by Colket (1986). The reaction of C 4 H 4 with the C 2 H 
radical was reported by Frenklach et al. (1983). Abstraction 
reactions of the oxygen atom and the hydroxyl radical with 
vinyl acetylene and those of C 2 H, the vinyl radical, and C 4 H 5 
with benzene were also used with estimated rate coefficients. 
Propargyl is removed by its direct reaction with molecular 


oxygen whose rate was measured by Slagel and Gutman 
(1986). Several of these reactions form 1,3 butadiene, which is 
removed by three reactions reported by Vaughn, Howard, and 
Longwell (1991). Recent discussions of benzene formation in 
the flames of aliphatic fuels (Miller and Melius, 1992 and 
Westmoreland, 1989) as well as a study of allene pyrolysis (Wu 
and Kern, 1987) have suggested the recombination of two 
propargyl radicals and the reaction of propargyl with allene as 
possible paths for the formation of benzene in aliphatic fuel 
oxidation. Several calculations with the rate coefficients given 
by Wu and Kern showed that these reactions had absolutely no 
effect on the computed results for benzene oxidation because 
propargyl is consumed very rapidly by its reaction with oxygen 
and its recombination with the H-atom to form allene. The two 
benzene-forming reactions were not included in the NASA 
mechanism. However, two additional reactions were added, 
both of which were used by Vaughn, Howard, and Longwell 
(1991) in their modeling of the benzene-ethylene-mixture 
oxidation in a perfectly stirred reactor: 


c 6 h 6 + ch 3 = c 6 h 5 + ch 4 

(10) 

C 6 H s + C 2 H 2 = C 6 H 5 C 2 H + H 

(11) 


Significant formation of the product phenylacetylene was 
observed and measured by Vaughn, Howard, and Longwell, 
who used only this one reaction to model their results. In the 
present work, reaction (11) formed phenylacetylene near the 
initial reaction temperature and then destroyed it as the tem- 
perature increased toward its equilibrium combustion value. 

One more point should be mentioned concerning reac- 
tion (20), the C 5 H 5 abstraction of a hydrogen atom from phenol 
(C 6 H 5 OH): 

C 5 H 5 + C 6 H 5 OH = C 6 H 5 0 + C 5 H 6 (20) 

The rate coefficient used in a previous work (Bittker, 1991) and 
in the present calculations is an estimate reported by Lovell, 
Brezinsky, and Glassman (1989) in their study of the pyrolysis 
mechanism of phenol; namely, 

k 2 Q = 2.67xl0 14 exp(-25200//?7) cm 3 /(mole-s) 

where R is the universal gas constant and T is the temperature in 
kelvin. The pre-exponential factor in this expression is rather high 
for this abstraction process and was replaced by another estimated 
expression in later work from the same research group. The latter 
rate coefficient was calculated from an estimate of the reverse 
reaction rate by means of the equilibrium constant and the law of 
microscopic reversibility. The expression, reported in the Emdee, 
Brezinsky, and Glassman (1992) paper is 

k 2Q = 4.20X10 13 r-°- 82 exp(-19840/#7) cm 3 /(mole-s) 
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This expression gives rate coefficient values about 2 orders of 
magnitude smaller than the first k 2Q given above by Lovell in 
1989. However, when the new expression was used in the 
behzene oxidation calculations shown in the next section, 
insignificant differences were found from all results computed 
using the older expression. These differences indicate that rate 
coefficients computed for reaction (20) by either expression 
are so high that the rate of this reaction cannot have any control 
over the initial stage of the benzene oxidation. There may be, 
however, situations in which reaction (20) becomes more rate 
controlling. Therefore, the newer and more realistic estimated 
k 2 Q is preferred and is listed in table II, which also includes the 
CH 4 , C 2 hydrocarbon, and H 2 oxidation reactions used in the 
original NASA Lewis mechanism. All reactions are allowed to 
be reversible, but only the forward rate coefficient parameters 
are given in the table because the reverse rates are generated by 
LSENS at every step of the reaction as the ratio of the forward 
rate coefficient to the equilibrium constant. This ratio was 
computed internally by LSENS from Gibbs function data as 
described previously. 

The NASA Lewis model differs from the Emdee model 
primarily in the products written for the reactions of molecular 
oxygen with the species benzene (C 6 H 6 ) and cyclopentadiene 

(C 5 H 6 ): 

C 6 H 6 + 0 2 = 0 6 H 5 O + OH (1) 

C 5 H 6 + 0 2 = C 5 H 5 0 + OH (22) 

In the Emdee mechanism, these reactions are simple 
hydrogen atom abstractions forming H0 2 and either the C 6 H 5 
or the C 5 H 5 radical. The above reactions, each of which could 
represent a two-step or more process, are an important source 
of OH radicals that initiate the reaction in the NASA Lewis 
mechanism. As discussed later in the next section, these differ- 
ences in products for reactions (1) and (22) are reflected in the 
results of sensitivity analysis computations to determine the 
reactions whose rates are most controlling in the entire oxida- 
tion process. 

Computational Results and Comparisons 

Comparisons of computed with experimental results were 
made with the expanded mechanism of 146 reactions. The 
experimental data are the same as those used by Bittker (1991). 
Discussed first are the concentration profiles measured by 
Lovell, Brezinsky, and Glassman (1988) in a highly mixed flow 
reactor at 101 kPa (1 atm) and at about 1 100 K. The benzene 
fuel was injected into a stream of nitrogen containing a low 
oxygen concentration and was mixed rapidly with the oxidant 
stream. The initial concentration and temperature conditions 
for the three cases reported are listed in table in as cases B-l, 
B-3, and B-4. 


Figures 1 and 2 show the computed concentration time 
profiles of the four species: benzene, carbon monoxide, 
cyclopentadiene, and phenol from the original NASA Lewis 
mechanism, the new NASA mechanism, and the experimental 
results of Lovell, Brezinsky, and Glassman (1988). The 
experiments were performed at three fuel equivalence ratios (p: 
0.14, 1.0, and 1.36. The time scale for the experimental 
concentration profiles was obtained by those authors from 
flow velocity measurements as a function of distance along the 
reactor and was corrected for the mixing time from the point of 
fuel injection into the oxidant stream. The flow reactor was 
modeled as a constant-pressure (101-kPa or 1-atm), homoge- 
neous, static reaction. Figure 1 shows that the new mechanism 
changes the computed results to reduce the overall reaction 
speed, as indicated by the time rate of benzene consumption 
and carbon monoxide formation. At (p = 0.74, the new com- 
puted results deviated more from experiment than did the 
original ones, which already showed a slower reaction than 
observed experimentally. However, at the other two equiva- 
lence ratios of 1.0 and 1.36, the new results are closer to the 
experimental curves because the original mechanism com- 
puted faster reaction rates, as measured by the benzene and 
carbon monoxide profiles. In figure 2, which shows phenol 
and cyclopentadiene composition versus time, the important 
observation is that the new mechanism significantly improved 
the prediction of the phenol and cyclopentadiene profiles at all 
three equivalence ratios. The concentrations of the radicals 
C 6 H 5 0 and C 5 H 5 are included in the computed concentrations 
of phenol and cyclopentadiene, respectively, because, as 
explained by Emdee, Brezinsky, and Glassman (1992), these 
resonantly stable radicals can achieve significant concentra- 
tions and could combine with a hydrogen atom in the sampling 
probe and be detected as the parent species. The original 
NASA Lewis mechanism predicted significantly more phenol 
formation than was observed for all experimental conditions. 
With the new mechanism, agreement is very good at all three 
equivalence ratios. The cyclopentadiene concentration pro- 
files computed for all conditions by the new mechanism are 
also in better agreement with the experimental results than 
were the original predictions. At the start of the reaction, the 
rate of increase of the C 5 H 6 concentration is still faster than the 
experimental rate is, but the computed maximum in this 
species concentration is significantly reduced. Compared with 
the original one, the new mechanism provides a marked 
improvement in the prediction of experimental results, except 
for the benzene and CO profiles at the lean equivalence ratio. 
This lack of improvement for these two species is most likely 
because there are still unknown reactions which make both 
mechanisms underpredict the benzene destruction rate (and 
corresponding carbon monoxide formation rate) at (p = 0.74 
and overpredict these rates at the other equivalence ratios. 

Experimental concentration profiles of six species are given 
in Emdee, Brezinsky, and Glassman (1992) for an additional 
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benzeneoxidation case in the same flow reactor. The equivalence 
ratio for this case is 0.91, and its initial conditions are listed as 
case B-2 in table III. Figure 3 shows the experimental benzene 
and phenol profiles for this case. Also in this figure are the 
computed profiles of Emdee’s work as well as computed 
profiles using the new mechanism of the present work. For 
benzene, the NASA Lewis computed profile is consistent with 
the experimental one but is slightly lower. The computed 
phenol concentration is in good agreement with experiment. A 
comparison of the NASA Lewis computed results with those of 
Emdee shows that their mechanism computes an excellent 
match to the observed benzene results, but the NASA Lewis 
mechanism computes a slightly more accurate phenol profile. 
In figure 4 are shown experimental and computed profiles using 
both mechanisms for carbon monoxide and C 2 hydrocarbons, 
which are mostly acetylene and a small amount of ethylene. The 
agreement between the NASA Lewis computation and the 
experimental results is good, but the computation gives higher 
concentrations than the experiment at the start of the reaction. 
However, the computed and experimental concentrations are 
quite close at the later reaction times shown. The Emdee 
mechanism overestimates the C 2 formation and underestimates 
the carbon monoxide formation. Experimental and computed 
concentration profiles of cyclopentadiene and C 4 hydrocar- 
bons (vinyl acetylene and 1,3 butadiene) are shown in figure 5. 
The NASA Lewis mechanism gives only a fair match to the 
experimental C 4 profile whereas Emdee’s mechanism gives 
numerically close results. The trend of the Emdee curve is, 
however, not consistent with that of the experimental profile. 
Neither Emdee’s nor the NASA Lewis benzene mechanism 
reproduces the cyclopentadiene concentration experimental 
profile. The NASA Lewis mechanism gives concentrations of 
C 5 H 5 + C 5 H 6 that are too high at early reaction times, and 
Emdee’s computations overestimate them for the entire reac- 
tion time interval shown in the figure. Both computed curves 
are, however, consistent with the experimental curve at long 
reaction times. These comparisons illustrate the strengths and 
weaknesses of each mechanism. 

The NASA Lewis mechanism was next tested on the predic- 
tion of the experimentally measured ignition delay times re- 
ported by Burcat, Snyder, and Brabbs (1986). Ignition was by 
reflected shock wave and the delay time was determined as the 
time (after shock passage) for observing the first significant 
pressure rise. For the computation, a static reaction at constant 
volume was assumed behind the shock, and the ignition delay 
time was taken as the time for a 5-percent pressure rise. Four 
benzene-oxygen-argon mixtures at three fuel equivalence 
ratios of 0.5, 1 .0, and 2.0 were used in the experiments. For the 
stoichiometric equivalence ratio, experiments covered two 
conditions, a dilute (95.616 mole % argon) and a strong 
(85.635 mole % argon) mixture. Initial conditions for all 
experiments are given in table IV along with the experimental 
results. Table IV also lists computed ignition delay times using 
the new mechanism and tabulates percentage differences and 


percent standard deviations <7 for each mixture. Ignition delay 
times t computed with the new mechanism are compared with 
the experimental results in figures 6 to 9 for all mixtures. These 
figures plot log 10 T versus the reciprocal of kelvin temperature 
T. The lines represent least-squares curve fits to the equation 

log 1Q r= Constant + AE/(2.3025S5*RT) 

where A E is an activation energy factor which is a measure of 
the temperature dependence of t, and R is the universal gas 
constant. Comparisons show that the present mechanism pre- 
dicts ignition delay times quite close to the experimental ones 
for the dilute mixture {(p = 1 .0) and gives fair agreement for the 
strong mixtures {(p- 0.5 and 1 .0). Agreement between experi- 
ment and calculation is poor for the (p = 2.0 mixture. Predicted 
ignition delay times are too long for the lean mixture and too 
short for both the strong stoichiometric and the rich mixtures. 
Similar results were shown in Bittker (1991) for the original 
benzene oxidation mechanism. This lack of agreement is attrib- 
uted to inadequacies in both mechanisms causing them to 
predict too slow a reaction for (p = 0.5 and too fast a reaction for 
the strong stoichiometric and rich mixtures. 

In view of the difficulty in defining the experimental first 
significant pressure rise, any <7 value less than 30 percent 
should be considered reasonable agreement between computa- 
tion and experiment. However, the present mechanism com- 
putes a weaker temperature dependence of ignition delay time 
than that observed experimentally for all mixtures. This fact 
can be seen from the slopes of the curve-fitted lines in fig- 
ures 6 to 9 and also from the activation energy factors A E 
computed from these slopes and shown in table V . The com- 
puted A E values are between 12 and 20 percent lower than the 
A E values for the experimental curve-fitted lines. 

In summary, the expanded benzene oxidation mechanism 
presented here has a significantly better capability of matching 
the available experimental concentration profile data in a flow 
reactor than does the original NASA Lewis mechanism. This 
is especially true for predicting phenol and cyclopentadiene 
concentration profiles at all experimental conditions. How- 
ever, the new mechanism is still incomplete and has some of 
the same inadequacies as the original mechanism in predicting 
higher temperature ignition delay times. The observed differ- 
ences between all experimental and predicted results will 
probably not be resolved until additional experimental results 
are available. 

Sensitivity Analysis Results 

The new benzene oxidation mechanism was developed with 
the help of detailed sensitivity analysis computations which 
determined the reactions whose rate coefficients had the great- 
est effect on computed results. Several reactions with a large 
uncertainty in their rate coefficients had the greatest effect on 
computed results when their coefficients were changed by a 
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moderate amount. These rate coefficients were adjusted to give 
the best overall agreement between computation and experi- 
ment, as shown in the previous section. Figure 10 presents 
normalized sensitivity coefficients for the reactions which 
control three concentration profiles shown in the previous 
section for the flow reactor case (p = 0.74. A normalized 
coefficient represents the approximate percent change in the 
variable that would be caused by a 1 -percent increase in the rate 
coefficient of a given reaction. A negative coefficient means 
that increasing the rate coefficient results in a decrease in the 
value of the dependent variable at that reaction time. These are 
local sensitivity coefficients and can only give information 
about moderate changes in the rate coefficient. They tell 
nothing about the effect of large changes in the latter and, 
therefore, cannot predict the effect of eliminating a reaction. 
The information in this figure will be compared with sensitivity 
analysis results for the Emdee mechanism. Emdee, Brezinsky, 
and Glassman (1992) list the following reactions as those 
whose rates have the strongest effect on the benzene profile: 


H + 0 2 = OH + O 

(130) 

c 6 h 5 o= c 5 h 5 + co 

(13) 

C 6 H 5 OH = C 6 H s O + H 

(17) 


The rate coefficient of their direct benzene-plus-oxygen 
reaction (a simple H-atom abstraction) has no effect on their 
computed results. The situation is quite different with the 
present NASA Lewis mechanism. 

Sensitivity coefficients for benzene and carbon monoxide 
concentrations are given in figures 10 (a) and (b). Reaction (1), 
the (chain branching) benzene reaction with molecular oxygen, 
is the major rate-controlling step of benzene consumption and 
corresponding formation of carbon monoxide. This situation is 
in contrast to the situation for the Emdee mechanism just noted. 
Reactions (13) and (17) are rate controlling in both mecha- 
nisms. The phenoxy dissociation (reaction (13)) speeds up the 
overall oxidation process whereas reaction (17), which pro- 
ceeds in reverse to form phenol, inhibits the oxidation. The 
H + 0 2 chain-branching step, which is the most rate-control- 
ling reaction in the Emdee mechanism, ranks much lower in the 
NASA Lewis mechanism, below two other rate-controlling 
steps, reactions (6) and (22), the OH attack on benzene and the 
reaction of cyclopentadiene with molecular oxygen. It must be 
emphasized that the H + 0 2 reaction is a very important reac- 
tion in the current mechanism, even though it is not the most 
rate-controlling one; because its rate is very rapid, it has a major 
effect on the radical pool concentrations, as it does in any 
hydrocarbon oxidation mechanism. 

Figure 10(c) shows sensitivity coefficients for the reactions 
which control phenol concentration. The main rate-controlling 
reactions for this species are two hydroxyl radical reactions, 
(6) and (23). Reaction (6) increases phenol formation whereas 


reaction (23), the destruction of phenol by OH, decreases 
phenol concentration. Net rate computations show that the 
reverse of reaction (17) is the dominant reaction that forms 
phenol, but this reaction ranks significantly below the last 
reaction in figure 10 for controlling the phenol concentration. 
Reaction (1) is fourth in the sensitivity ranking for this species. 

To find out what reactions are rate controlling on the ignition 
delay times shown in the previous section, sensitivity coeffi- 
cients of pressure were computed and are shown in figure 1 1 for 
the shock ignition of the strong stoichiometric mixture at 
1435 K. The same three reactions which control benzene 
destruction and carbon monoxide formation in the lower tem- 
perature flow reactor experiments, (1), (13), and (17), are seen 
to control the pressure rise in the reflected shock ignition 
reactions. As would be expected, they also control benzene and 
carbon monoxide concentration profiles computed for this 
ignition reaction. The H + 0 2 chain-branching reaction is next 
in importance (at the higher temperature of this reaction) for 
controlling pressure and also carbon monoxide concentration 
(sensitivity coefficients not shown). This is a significant increase 
in its rate-controlling effect compared with its sensitivity rank- 
ing for the lower temperature flow reactor experiments. 

Toluene Oxidation 

Chemical Mechanism 

The toluene oxidation model presented by Emdee, Brezinsky, 
and Glassman (1992) was added to the new benzene model just 
discussed. Included are 41 of the first 42 reactions in their 
table IV, which are reactions of toluene and its fragments. The 
one reaction omitted is the slower of two paths given for the 
reaction of the benzyl radical (C 6 H 5 CH 2 ) with the H0 2 radi- 
cal. Computations with and without the slower process gave 
identical results, so it was not used. From their table, four 
additional reactions involving the species formaldehyde and 
CH 2 OH were also used. Other reactions in their table were 
included in the new NASA benzene model. The 45 reac- 
tions in the toluene model and their rate coefficient param- 
eters are listed in table VI. Preliminary calculations were 
made using all the Emdee rate coefficient parameters to 
compute the experimental concentration profiles for the 
two toluene flow reactor cases shown in the Emdee paper. 
The experimental conditions of Burcat, Snyder, and Brabbs 
(1986) were used to also compute the pressure-defined igni- 
tion delay times for toluene oxidation behind a reflected shock. 
Sensitivity coefficients were also computed to determine which, 
if any, of the Emdee rate coefficients might need adjusting to 
obtain the best agreement between computed and experimen- 
tal results. 

The preliminary sensitivity analysis showed that only five 
of the Emdee rate coefficients had to be adjusted to achieve the 
best overall agreement between the computations and all the 
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experimental results. In agreement with Emdee, the main rate- 
controlling reaction was the direct toluene-oxygen reaction 

C 6 H 5 CH 3 + 0 2 = C 6 H 5 CH 2 + H0 2 (147 ) 

Because the rate coefficient of this reaction has not been 
measured experimentally, its pre-exponential factor was used 
as one of the important adjustable parameters for matching the 
experimental results. The value of 2.5 X 10 14 given in table VI 
for A 147 is only 16.7 percent lower than the Emdee value 
obtained by this same procedure. Two other reactions whose 
pre-exponential factors had to be changed from Emdee ’ s values 
are the reactions of the H-atom with toluene: 

C 6 H 5 CH 3 + H = C 6 H 6 +CH 3 (152) 

C 6 H 5 CH 3 + H = C 6 H 5 CH 2 + H 2 (153) 

The rates of both these reactions have a significant effect on the 
concentration time profiles of benzene, methane, and styrene. 
These three species profiles were matched to the experimental 
results by increasing A 152 to 25 percent above Emdee’s value 
and decreasing A 153 to 16.7 percent below Emdee’s value. 
Reaction (153) also influences the cresol profile, and the 
decrease in its A-factor also improved this species’ agreement 
between computed and experimental results. Only two other 
pre-exponential factors had to be changed from Emdee’s 
values: these are the A-factors for the reactions 

CH 3 C 6 H 4 0 = C 6 H 6 + H + CO (157) 

and 

C 6 H 5 CH 2 + O = C 6 H 5 CHO + H (164) 

Reaction (157) controls the cresoxy concentration, which con- 
trols the formation of cresols by the reaction 

CH 3 C 6 H 4 0 + H = CH 3 C 6 H 4 OH (156) 

The value of A 157 had to be decreased by 74 percent from 
Emdee’s value to increase the computed cresol concentration 
level to values near those reported in the experiment. This 
change was justified, inasmuch as Emdee’s rate coefficient 
expression is not an experimentally measured one. It is the 
expression determined experimentally by Lin and Lin (1986) 
for the analogous decomposition of thephenoxy radical C 6 H 5 0 
(see reaction (13) in table II). Reaction (164) is an important 
rate-controlling reaction for benzaldehyde. The estimated pre- 
exponential factor of Emdee was increased by 50 percent to 
give the best agreement between computed and experimental 
benzaldehyde concentration profiles. 

In the next section are shown comparisons of the new NASA 
computed results with the experimental results of Emdee, 


Brezinsky, and Glassman and Burcat, Snyder, and Brabbs for 
toluene oxidation using this modified Emdee toluene mecha- 
nism coupled with the new benzene oxidation mechanism. 

Computational Results and Comparisons 

As indicated earlier, two sets of experimental data for high- 
temperature toluene oxidation were used in the development of 
the NASA complete oxidation mechanism. The first set con- 
sists of concentration profiles measured in the same highly 
mixed, atmospheric-pressure flow reactor used to obtain the 
benzene oxidation data described earlier in this paper. The two 
cases reported by Emdee, Brezinsky, and Glassman (1992) had 
the initial conditions given in table III as cases T-l and T-2. Of 
course, the computer model of these experiments is the same as 
that for the benzene experiments. The second set of data 
consists of ignition delay times for toluene-oxygen-argon- 
mixture ignition behind a reflected shock wave. These data are 
part of the same work of Burcat, Snyder, and Brabbs (1986), 
who reported the benzene-oxygen ignition delay times 
described earlier in the section Benzene Oxidation. The 
pressure-rise-based ignition delay times were calculated 
exactly as were the ones for benzene. The flow reactor experi- 
ments are discussed first. 

Figures 12 (a) and (b) present experimental concentration 
profiles of benzene and toluene for the lean {(p- 0.69) and rich 
((p= 1.33) cases of table III, respectively. Also plotted in this 
figure and in figures 1 3 to 1 8 are two computed curves for each 
concentration profile. These are the curves using the Emdee 
mechanism (from their paper) and the curves using the oxida- 
tion mechanism of the present work. Figure 1 2 shows that both 
the Emdee and NASA Lewis mechanisms give excellent agree- 
ment with experiment , except for the benzene profile in the lean 
case for which the NASA Lewis curve is slightly low. Concen- 
tration profiles for benzaldehyde (C 6 H 5 CHO) and methane are 
presented in figure 13 for the case (p = 0.69. For each of these 
species, both mechanisms give good agreement with the 
experimentally observed trends and give fairly good quantita- 
tive agreement. The computed benzaldehyde profiles slightly 
over predict the experimental concentrations for early reaction 
times but then come into good agreement with the observed 
profile. The computed maximum concentration and corre- 
sponding reaction time for both computed curves are in good 
agreement with the experimental values, and the slopes of the 
computed and experimental curves are in good agreement. The 
NASA Lewis mechanism reproduces the experimental meth- 
ane profile quite well and somewhat better than the Emdee 
computation, which underpredicts its concentration. Figure 14 
presents experimental and computed benzyl alcohol 
(C 6 H 5 CH 2 OH) concentrations for the same case. The agree- 
ment between computed and experimental profiles is slightly 
better for the Emdee mechanism. For both computed curves, 
the maximum concentrations are in good agreement with 
experiment, but the times for reaching this concentration are 
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too long and the observed sharp drop in concentration from the 
peak value is not well predicted. Figure 15 plots experimental 
and computed benzyl alcohol and benzaldehyde concentra- 
tions for the case (p - 1 .33. Both mechanisms do a reasonably 
good job of predicting both profiles. The matching is better for 
benzyl alcohol than for benzaldehyde and the Emdee mecha- 
nism gives a more accurate maximum benzyl alcohol 
concentration. 

Comparisons of computed and experimental concentration 
profiles for other species are given in figures 16 to 1 8. They are 
all for the case <p= 0.69. Figure 16 shows excellent agreement 
between computed and experimental carbon monoxide profiles 
for the NASA computed curve. Emdee’ s computed concentra- 
tion curve for this species also matches the experimental curve 
quite well but shows some underprediction of the experimental 
results at reaction times greater than 70 ms. The acetylene 
concentration profiles in figure 16 show that the new mecha- 
nism consistently overpredicts acetylene concentration. At 
reaction times greater than 75 ms, the computed concentration 
is 1.5 to 2 times higher than the experimental value. By 
comparison, Emdee’ s computed acetylene profile also 
overpredicts the experimental results and shows much more 
deviation than the NASA computed curve. At reaction times 
greater than 75 ms, the Emdee computed acetylene concentra- 
tions are approximately three times the experimental values. 
Figure 17 shows that concentration profiles for phenol, 
ethylbenzene, and styrene are fairly well predicted by the 
NASA mechanism. The observed peak in phenol concentration 
is, however, not well predicted by the present computations. 
The Emdee computations predict ethylbenzene and styrene 
concentrations about as well as the NASA mechanism does. 
However, Emdee’ s computed phenol concentrations are sig- 
nificantly below the level of the experimental curve. As was 
done for the benzene oxidation computations, the experimental 
phenol concentration was compared with the sum of the com- 
puted phenol and phenoxy radical concentrations. The phenoxy 
concentration makes only a small contribution to the total of 
both concentrations. The poorest matching of experimental and 
computed concentration profiles is shown in figure 18 for the 
species cyclopentadiene and cresol. Again, the concentrations 
of the cresoxy and cyclopentadienyl radicals were included in 
the cresol and cyclopentadiene concentrations. Both are sig- 
nificantly underpredicted by the NASA Lewis mechanism. 
However, the figure shows that the NASA computed cresols 
curve goes through a maximum concentration value at reaction 
time of 120 ms, which is quite close to the experimental 
concentration maximum at 75 ms. Also, after a reaction time of 
40 ms, the slope of the NASA computed cresol curve is about 
the same as that of the experimental curve. A shift in the time 
scale of the experimental curve would bring it into much better 
agreement with the computed curve. The Emdee mechanism 
underpredicts cresols noticeably more than the NASA Lewis 
mechanism does. It is interesting to note, however, that the 
Emdee curve for C 5 H 5 + C 5 H 6 overpredicts the measured 


cyclopentadiene concentration. Their computed concentration 
profile for C 5 H 6 alone (not shown) matches the experimental 
curve fairly well. 

The NASA toluene mechanism was next used to compute the 
pressure-rise-based ignition delay times for the toluene-oxygen- 
argon mixtures mentioned at the beginning of the Toluene 
Oxidation section. Both the experimental values of Burcat, 
Snyder, and Brabbs (1986) and the NASA computed values 
were obtained exactly as described previously for the shock- 
initiated ignitions of benzene and oxygen. Four different mix- 
tures and two equivalence ratios (p= 0.33 1 and 1 .0) were used 
by Burcat. The mixture descriptions are given in table VII, 
which also tabulates the experimental and computed delay 
times rand gives the same error analysis shown in table IV for 
the benzene ignitions. Computed and experimental ignition 
delay times are plotted as the log 10 f versus the reciprocal of 
temperature for each mixture in figures 19 to 22. In each figure 
are plotted the actual experimental data points and the least- 
squares line fitted through the points. For the computations, 
only the least-squares-fitted lines through the computed delay 
times are shown. Good agreement between experimental and 
computed results is shown in figures 20 and 21 for mixtures 2 
and 3. Both these mixtures are stoichiometric and were diluted 
with 95.027 mole % argon. Mixture 2 experiments had initial 
pressures of approximately 202 kPa (2 atm) whereas mix- 
ture 3 cases had initial pressures close to 606 kPa (6 atm). 
Table VII shows that the standard deviations of the computed 
rvalues are 10.7 and 13.1 percent, respectively, for mixtures 2 
and 3. Figure 19 shows that the computed rvalues are sig- 
nificantly longer than the experimental ones for mixture 1 
((p = 0.331; initial pressures ~ 202 kPa or 2 atm). In contrast, 
computed r values are much shorter than the experimental 
ones for mixture 4, which is stoichiometric but diluted 
in only 85.053-mole % percent argon. These results are shown 
in figure 22. Table VII gives the standard deviations of 42. 1 and 
36.4 percent for mixtures 1 and 4, respectively. The poor 
agreement for mixture 1 can be explained by the possibility of 
as yet undiscovered reactions which are important at the lean 
equivalence ratio and not at (p = 1 .0. However, this explanation 
does not apply to the poor agreement for mixture 4. The only 
difference between mixtures 2 and 4 is the amount of argon 
dilution. The enthalpies of the mixture 4 cases at the unshocked 
temperature of approximately 298 K are significantly higher 
than those of the mixture 2 cases because of the larger fraction 
of fuel and oxidant in mixture 4, which also has higher heat 
capacity. Therefore, one would expect the initial temperature 
range of the mixture 4 cases to be lower than that of the mixture 
2 cases for the same range of reflected shock pressure. T able VII 
shows that this is true. However, it is possible that these 
temperatures, which were computed using the idealized shock 
equations, are slightly high for mixture 4. A small estimated 
constant correction of about 0.2 percent for shock velocity 
attenuation was applied by this author to Burcat’ s original 
calculations of reflected shock conditions. If this correction 
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were greater for the higher enthalpy mixture or other deviations 
from ideal gas conditions became important, computed tem- 
peratures could be lowered. A modest lowering of even 
1 percent (14 to 15 K) produces significantly improved agree- 
ment between computed and experimental results for the mix- 
ture 4 ignition delay times. The standard deviation of all points 
is reduced from 36.4 to 26.6 percent, with maximum and 
minimum deviations of about 34 and 15 percent for the indi- 
vidual data points. Also to be mentioned, as done for the 
benzene ignitions, is the uncertainty in determining the time of 
the first significant pressure rise from the experimental pres- 
sure time traces. Considering these experimental uncertainties, 
this comparison of computed and experimental ignition delay 
times raises questions that will be answered by improving the 
chemical mechanism and by obtaining more definitive experi- 
mental measurements of ignition delay times and other data for 
the oxidation of toluene at temperatures above 1300 K. 

A comparison of the NASA computed temperature depen- 
dence of ignition delay and the experimental observations is 
shown in table VIII, which lists the activation energy factors A E 
defined for the benzene ignition delay measurements. They 
were computed from the slopes of the least-squares lines shown 
in figures 19 to 22. As in the case of benzene, the general trend 
is that the computed A E values are lower than the experimental 
values, in this case by about 10 to 19 percent. However, for 
these toluene mixtures, only three of the four are included in this 
trend whereas all mixtures for benzene followed the trend. 
Table VIII shows that the computed A E for the dilute and high- 
pressure stoichiometric mixture 3 is about 23 percent higher 
than the experimental value. This result is attributed to the low 
experimental A E, which is about 37 percent lower than the A E 
values for mixtures 1 and 2. Certainly, one would expect the 
temperature dependences of the three stoichiometric mixtures 
to be close to each other and also be close to that of the lean 
mixture. Table V shows that all benzene mixtures had about the 
same experimental temperature dependence. This discrepancy 
also shows the need for additional higher temperature data in 
ignition delay times for toluene-oxygen mixtures. 

In summary, the NASA Lewis combined toluene-benzene 
oxidation mechanism presented herein accurately predicts most 
of the concentration profiles measured for experiments in a 
highly mixed flow reactor at about 1200 K. When it is used to 
predict ignition delay times for higher temperature (1300 to 
1600 K) shock-initiated oxidations, the results are mixed. The 
computed pressure-based ignition delay times agree well with 
experimental results for some conditions but give only poor to 
fair agreement for other conditions. As was the case with the 
benzene oxidation mechanism alone, as yet undiscovered reac- 
tions must be added to the mechanism and also more experi- 
mental results are needed to resolve the discrepancies observed 
between the presently available experimental data and com- 
puted results, especially at temperatures above 1300 K. It 
should be mentioned that the Emdee mechanism also failed to 
give good agreement with experiment when investigators applied 


it to the calculation of laminar flame speeds for both benzene 
and toluene mixtures. Computed flame speeds for this higher 
temperature regime were lower than reported experimental 
values (Emdee, Brezinsky, and Glassman, 1992). 

Sensitivity Analysis Results 

Some typical results of the NASA Lewis sensitivity analysis 
computations performed during the development of the toluene 
oxidation mechanism are shown in figure 23. Sensitivity coef- 
ficients of reactions that control four species profiles are shown 
for case T- 1 of table III, the lean-mixture toluene oxidation. The 
coefficients were all computed at a reaction time of 60 ms. 
These results are in general agreement with those of Emdee, 
Brezinsky, and Glassman (1992) for the reactions which con- 
trol the toluene concentration profile (fig. 23(a)). The most 
sensitive steps are reactions (147) and (130), which have 
negative sensitivity coefficients and therefore promote the 


consumption of toluene: 



c 6 h 5 ch 3 + o 2 

= c 6 h 5 ch 2 + ho 2 

(147) 

H + 0 2 

= OH + O 

(130) 


It should be mentioned that reaction (147) consumes only a 
small fraction of the toluene. It is a very sensitive step because 
it is the main source of the H0 2 radical, which is then consumed 
by reaction with the benzyl radical to form benzaldehyde, H, 
and OH radicals (reaction (176)). Net species formation rates 
show that reaction (176) is a primary source of the H-atom, 
which propagates the chain reaction destruction of toluene. 
Figure 23(a) also shows that H0 2 reacts with the HCO radical 
in reaction (118), the sixth most sensitive reaction with a 
negative sensitivity coefficient for toluene. Its products are 
molecular oxygen and formaldehyde, CH 2 0, which pyrolyzes 
to form H-atoms and more HCO. Both products, therefore, 
promote the chain reaction oxidation of toluene. The next two 
most sensitive rate coefficients are for the reactions 

C 6 H 5 CH 3 + H = C 6 H 5 CH 2 + H 2 (153) 

and 

H + C 6 H 5 CH 2 = C 6 H 5 CH 3 (148) 

Both of these reactions proceed from left to right as written. 
Although reaction (153) does destroy toluene, it also forms the 
benzyl radical, which participates in reaction (148), the most 
important reaction that reforms toluene and would have an 
inhibiting effect on the oxidation. Figure 23(a) shows, how- 
ever, that the benzyl radical also reacts with benzaldehyde in 
the fifth most sensitive step, reaction (187). This reaction has a 
negative sensitivity coefficient for toluene and, thus, an accel- 
erating effect on the oxidation. Even though this reaction forms 
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toluene, the major effect of this reaction is the formation of 
C 6 H 5 CO, which rapidly dissociates to form CO and ultimately 
generates other free radicals that speed up fuel consumption 
and product formation. The net result of these competing 
qualitative effects in a complex mechanism usually has to be 
determined by a detailed sensitivity analysis. The positive 
sensitivity coefficient of reaction (153) in figure 23(a) shows 
that its overall effect is to inhibit the chain reaction process that 
consumes toluene. The dominance of the inhibition can be 
explained by the following observations. First, net-species- 
formation-rate calculations show that reaction (153) is the 
largest consumer of H-atoms. It consumes more H-atoms than 
the important chain-branching step, H + 0 2 , reaction (130), and 
so competes with it and slows down the chain oxidation of 
toluene. Second, the benzyl radical formed by reaction (153) 
speeds up reaction (148), which has two inhibiting influences: 
it not only consumes more H-atoms but also reforms toluene. 
The positive coefficient of reaction (148) is consistent with its 
primary effect of forming toluene. Therefore, both these reac- 
tions inhibit the oxidation process if their rate coefficient is 
increased. The Emdee sensitivity analysis lists the reaction 

C 6 H 5 CH 3 + OH = C 6 H 5 CH 2 +H 2 0 (154) 

as the fifth most sensitive reaction for toluene concentration 
and its sensitivity coefficient is positive. The present analysis 
also found a positive sensitivity coefficient for this reaction, but 
it ranks 10th in sensitivity, well below the last four reactions 
listed in figure 23(a). In the oxidation of aliphatic hydrocar- 
bons, radical attacks on the fuel molecule (like reactions (153) 
and (154)) usually increase fuel consumption. However, both 
these reactions inhibit toluene consumption because of the 
further reaction of the benzyl radical just discussed. 

The reactions that control toluene also control CO concen- 
tration, as shown by the sensitivity coefficients in 
figure 23(b). Reactions (130) and (147) are the most important 
in promoting CO formation, followed by reaction ( 1 87). Reac- 
tions (148) and (153) are the primary inhibtitors of CO forma- 
tion. It should be noted that two benzene oxidation reactions are 
also in the group of the eight most sensitive reactions control- 
ling both toluene and CO concentration. The first is the phenoxy 
dissociation (reaction (13)) which helps consume toluene and 
also form CO. The other is reaction (17), which reacts in the 
reverse direction to consume phenoxy and a hydrogen atom in 
forming phenol. This reaction, therefore, inhibits the oxidation 
process by competing for H-atoms with reaction (130), as 
shown by its positive sensitivity coefficient for toluene concen- 
tration and its negative coefficient for CO concentration. 

Reactions (130) and (147) are dominant rate-controlling 
reactions in the formation of cresols and benzyl alcohol, as 
shown in figures 23(c) and (d), respectively. They increase the 
formation of both species. In figure 23(c) the next two most 
sensitive reactions inhibit the formation of cresols. These are 
reactions (153) and (154), the H and OH attacks on toluene. 


Note that the OH reaction is much more rate controlling 
on cresols than on the CO and toluene concentrations. 
Reaction (187) is again rate controlling and increases cresol 
concentration. Reaction ( 1 57), the decomposition of the cresoxy 
radical, is the sixth-ranked rate-controlling step. This reaction 
has a negative sensitivity coefficient because removing cresoxy 
reduces the generation of cresols by its main path, the combi- 
nation of an H-atom with cresoxy (reaction (156)). These 
results of the sensitivity analysis led to the decreasing of the 
Emdee pre-exponential factor of reaction ( 1 57) to significantly 
increase cresol concentration without changing any other 
species concentrations. Although figure 23(c) shows that reac- 
tion (156) itself is the eighth-ranked controlling step, there was 
no attempt to increase its rate coefficient above the Emdee 
value. This value is already high and it would have had to have 
been increased by a very large amount to effect the cresol 
concentration significantly. This Emdee rate coefficient is 
essentially the value measured by He, Mallard, and Tsong 
(1988) for the analogous combination of the phenoxy radical 
with the H-atom to form phenol. 

From figure 23(d), one sees that the most important reaction 
in controlling benzyl alcohol concentration is the principal path 
for its formation, reaction (165); however, the three previous 
sensitivity coefficient plots indicate that this is not usually the 
case. The rate coefficient for this recombination of the OH 
radical with the benzyl radical was used by Emdee, Brezinsky, 
and Glassman as an adjustable parameter to match the com- 
puted and experimental benzyl alcohol concentration profiles; 
the same procedure was used in the present work. After the rate 
coefficient of reaction (147) was set to obtain the best agree- 
ment with all other concentration profiles, the coefficient of 
reaction (165) was adjusted to give the best benzyl alcohol 
profile without changing any of the other concentration pro- 
files. It is not surprising that the NASA Lewis fc 165 agreed 
exactly with Emdee’ s value, inasmuch as his thermodynamic 
data for the benzyl radical and the benzyl alcohol species were 
used. It was observed here that the attack of the OH radical on 
toluene, reaction (154), is the main inhibitor of benzyl alcohol 
formation. Reaction (153) the attack of the H-atom on the fuel, 
is a less sensitive reaction and is now a reaction which promotes 
the formation of benzyl alcohol. This reaction inhibits the 
formation of the three other species already discussed. The 
reason why reactions (153) and (154) have different effects on 
benzyl alcohol concentration is not clear because they both 
form the benzyl radical and thus help reaction (165) to form 
benzyl alcohol. One possibility is that the water formed in 
reaction (154) is much less reactive than the molecular hydro- 
gen formed in reaction (153). This example shows again that it 
is not always clear what the effect of changing a given reaction 
rate coefficient will be in a complex system, and a detailed 
sensitivity analysis needs to be done. 

To find out which reaction rates control the computed 
toluene oxidation ignition delay times presented in the previous 
section, the sensitivity coefficients of pressure were calculated 
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for the present toluene mechanism. Typical results are given in 
figure 24 for the case (p - 0.331 at an initial temperature of 
1334 K. The reaction time of 700 ms is about 200 ms before the 
ignition delay time. The same reactions which control the flow 
reactor oxidation experiments can be seen to control the shock 
ignitions at the higher temperatures (and pressures). The pres- 
sure rise is controlled mainly by reactions (130) and (147) and 
is inhibited by the H and OH attacks on the fuel, reactions (153) 
and (154). The dissociation-recombination reactions of phenol 
and phenoxy from the benzene oxidation mechanism also 
influence the rate of pressure rise. Although not shown, the 
sensitivity analysis results for toluene concentration give the 
same reaction rankings as those shown for pressure. 

In summary, the sensitivity analysis findings were in general 
agreement with the results of Emdee, Brezinsky , and Glassman 
(1992) but differed in some details. The two primary rate- 
controlling steps were the direct toluene-oxygen reaction and 
the H + 0 2 radical chain-branching step. The inhibiting toluene 
+ OH reaction was less controlling on the toluene and also CO 
concentration profiles than Emdee found. In the NASA Lewis 
mechanism, the reaction of benzaldehyde with benzyl, reaction 
(187), was a significant rate-controlling step. It promoted the 
consumption of fuel and the corresponding formation of CO. 
The present analysis also showed that two benzene mechanism 
reactions had a noticeable effect on the computed concentration 
and pressure rise profiles in toluene oxidation. 

Conclusions 

This report presented an expanded NASA Lewis benzene 
oxidation mechanism which has improved capability to match 
most experimental concentration profiles. This mechanism 
was compared with that developed by Emdee, Brezinsky, and 
Glassman and differs from it in several respects. The two 
mechanisms were comparable in computing published experi- 
mental concentration profiles. For example, the Emdee mecha- 
nism predicted more accurate benzene and C 4 hydrocarbon 
profiles whereas the NASA Lewis mechanism gave more 
accurate CO, phenol, and C 2 hydrocarbon profiles. Neither 
mechanism predicted the cyclopentadiene concentration pro- 
file very well. The NASA Lewis mechanism also did a reason- 
able job of computing higher temperature and pressure 
experimental ignition delay times. A modified version of the 
Emdee toluene oxidation mechanism was coupled to the NASA 
Lewis benzene mechanism. A few of the Emdee rate coeffi- 
cients had to be modified because of differences between the 
two benzene submodels. The NASA Lewis complete model 
was used to compute not only Emdee’ s toluene experimental 
data but also to match higher temperature and pressure toluene 
oxidation ignition delay times. Comparisons of computed 
results for both complete toluene models with experimental 
data showed that each was equally good in matching the 


experiments. Both had successes and failures in reproducing 
the experimental results. When the NASA Lewis predictions of 
toluene oxidation concentration profiles were compared with 
those of the Emdee mechanism, several species concentrations 
were predicted equally well. However, Emdee matched the 
experimental benzyl alcohol profiles better than the NASA 
Lewis mechanism did; the NASA Lewis matched the experi- 
mental methane profile better than Emdee did. Although both 
mechanisms predicted more acetylene production than observed 
experimentally, the NASA computed curve was noticeably 
closer to the experimental one. Also, the NASA Lewis com- 
puted phenol concentration profile agreed much better with the 
experimental results than did Emdee’ s, which was well below 
the experimental curve. Both mechanisms had difficulty com- 
puting the experimental cresol concentration profile. One ma- 
jor difference between the two mechanisms was their computed 
concentration profiles of cyclopentadiene. The NASA Lewis 
mechanism predicted too low concentrations of this species and 
the Emdee mechanism, too high levels. 

When both benzene and toluene oxidation mechanisms 
(Emdee and NASA Lewis) were applied to matching experi- 
mental data taken at temperatures above 1200 K, mixed results 
were obtained, i idicating that they both lacked important, as 
yet undetermined reactions needed to fully explain high- 
temperature experimental results. The Emdee mechanism com- 
puted laminar flame speeds which are lower than values 
reported in the literature for both benzene and toluene flames. 
In the present work, both benzene and toluene mechanisms 
gave only mediocre agreement when used to compute experi- 
mental ignition delay times for oxidations behind a reflected 
shock. These mechanisms must be applied to new high- 
temperature experimental results as they become available to 
obtain a complete oxidation mechanism for both fuels. It is also 
clear that several uncertain reaction rate coefficients for key 
steps in these mechanisms need to be measured. 

In conclusion, both the NASA Lewis and Emdee models for 
toluene and benzene oxidation are comparable in predicting the 
presently available experimental data for these fuels. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 31, 1995 
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Parameter 


Temperature, 

T, 

K 


Phenol 

C 6 H 5 OH 


TABLE I. — THERMODYNAMIC DATA FOR REACTING SPECIES 


Species 


Phenoxy Benzene 

c 6 h 5 o c 6 h 6 


(a) Thermodynamic data at 298.15 K 


Cyclopentadienolyl 

C 5 H 4 OH 


Cyclopentadienonyl 

C 5 H 5 0 


Heat capacity, C p , 
cal/mole-K 

24.789 

22.566 

19.629 

18.873 

22.130 

21.516 

Enthalpy, A jH, 
kcal/mole 

-23.034 

11.404 

19.792 

78.503 

19.150 

19.300 

Entropy, S°, 
cal/mole-K 

75.329 

73.570 

64.308 

68.914 

73.150 

73.541 


(b) Heat capacity C D and entropy S° at various temperatures, cal/mole-K 


400 

32.388 

83.705 

29.766 

81.234 

26.924 

71.121 

25.453 

75.403 

28.819 

80.617 

28.050 

600 

43.685 

99.134 

40.611 

95.501 

38.099 

84.293 

35.527 

87.756 

38.505 

94.285 

37.938 

1000 

55.571 

124.61 

52.169 

119.33 

50.419 

107.02 

46.590 

108.84 

48.285 

116.54 

48.052 

1200 

58.679 

135.02 

55.087 

129.10 

53.624 

116.51 

49.420 

117.60 

50.789 

125.57 

50.679 

1400 

61.261 

144.27 

57.497 

137.78 

56.278 

124.98 

51.761 

125.40 

52.820 

133.56 

52.823 

1600 

63.381 

152.59 

59.462 

145.59 

58.448 

132.64 

53.671 

132.44 

54.443 

140.72 

54.548 

1800 

65.097 

160.16 

61.041 

152.69 

60.197 

139.63 

55.208 

138.85 

55.717 

147.21 

55.916 

2000 

66.465 

167.09 ! 

62.288 

159.19 

61.585 

146.05 

56.424 

144.73 

56.699 

153.13 

56.981 


80.796 

94.182 

116.23 

125.23 
133.21 
140.38 
146.88 
152.83 























































TABLE I.— Continued. THERMODYNAMIC DATA FOR REACTING SPECIES 


Parameter 

Species 


Cyclopentadiene 

c 5 h 6 

Cyclopentadienyl 

c 5 h 5 

Butadiene 

c 4 h 6 

Butadienyl 

C 4 H 5 

Ketyl radical 

c 2 ho 

Ketene 

c 2 h 2 o 

(a) Thermodynamic data at 298.15 K 

Heat capacity, C p , 
cal/mole-K 

18.041 

18.309 

19.020 

18.317 

10.803 

12.367 

Enthalpy, AjH, 
kcal/mole 

31.935 

63.500 

26.109 

76.007 

38.500 

-11.401 

Entropy, S°, 
cal/mole-K 

65.495 

66.798 

66.630 

70.585 

58.956 

57.814 


(b) Heat capacity C p and entropy S° at various temperatures, cal/mole-K 


Temperature, 

T, 

K 


B 


■ 


■ 


■ 

C /> 

■ 


■ 

400 

24.646 

71.744 

24.686 

73.099 

24.606 

73.016 

23.224 

76.678 

11.716 

62.265 

14.215 

61.717 

600 

34.697 

83.768 

33.484 

84.925 

32.746 

84.663 

30.253 

87.529 

13.068 

67.282 

16.885 

68.020 

1000 

45.794 

104.42 

42.690 

104.48 

41.467 

103.70 

37.697 

104.97 

15.243 

74.490 

20.245 

77.509 

1200 

48.682 

113.03 

45.402 

112.51 

44.073 

111.49 

39.724 

112.03 

16.089 

77.347 

21.356 

81.303 

1400 

51.060 

120.72 

47.412 

119.67 

46.033 

118.45 

41.346 

118.28 

16.767 

79.880 

22.205 

84.662 

1600 

52.992 

127.67 

48.925 

126.10 

47.413 

124.69 

42.626 

123.89 

17.302 

82.155 

22.855 

87.671 

1800 

54.540 

134.01 

50.081 

131.94 

48.497 

130.34 

43.619 

128.97 

17.723 

84.218 

23.359 

90.393 

2000 

55.761 

139.82 

50.977 

137.26 

49.354 

135.50 

44.378 

133.60 

18.057 

86.104 

23.752 

92.875 



























































TABLE I.— Continued. THERMODYNAMIC DATA FOR REACTING SPECIES 


Parameter 

j Species j 


Toluene 

Benzyl 

Benzyl alcohol 

Cresols 

Ethyl benzene 

Styrene 


c 7 h 8 

c 6 h 5 ch 2 

C 6 H 5 CH 2 OH 

ch 3 c 6 h 4 oh 

c 6 h 5 c 2 h 5 

c 6 h 5 c 2 h 3 


(a) Thermodynamic data at 298.15 K 


Heat capacity, C p , 
cal/mole-K 

24.684 

26.082 

26.389 

30.599 

30.449 

28.726 

Enthalpy, AjH, 
kcal/mole 

11.991 

50.311 

-24.000 

-31.620 

7.151 

35.445 

Entropy, S°, 
cal/mole-K 

76.527 

76.753 

80.069 

86.072 

86.193 

82.402 


(b) Heat capacity C p and entropy S° at various temperatures, cal/mole-K 


Temperature, 

T, 

K 


■ 

D 

■ 


S° 


S° 


■ 

S 

■ 

400 

33.105 

84.970 

34.030 

85.550 

35.745 

89.144 

39.439 

96.326 

40.348 

96.535 

37.943 

92.193 

600 

46.927 

101.16 

46.468 

101.86 

51.061 

106.71 

53.077 

115.07 

56.774 

116.18 

52.521 

110.48 

1000 

62.452 

129.19 

60.224 

129.20 

67.559 

137.12 

67.964 

146.09 

75.117 

149.98 

68.269 

141.43 

1200 

66.668 

140.96 

63.928 

140.52 

71.928 

149.84 

72.016 

158.85 

80.155 

164.14 

72.433 

154.26 

1400 

70.114 

151.50 

66.954 

150.61 

75.500 

161.21 

75.336 

170.21 

84.230 

176.82 

75.822 

165.69 

1600 

72.895 

161.05 

69.393 

159.71 

78.382 

171.48 

78.021 

180.45 

87.522 

188.29 

78.543 

176.00 

1800 

75.105 

169.77 

71.331 

168.00 

80.673 

180.85 

80.162 

189.77 

90.146 

198.75 

80.695 

185.38 

2000 

76.833 

177.78 

72.845 

175.60 

82.466 

189.45 

81.843 

198.31 

92.206 

208.36 

82.369 

193.97 




















































TABLE I.— Concluded. THERMODYNAMIC DATA FOR REACTING SPECIES 


Parameter 

Species 


Cresoxy 

Benzaldehyde 

Hydroxyl 


ch 3 c 6 h 4 o 

c 6 h 5 cho 

OH 


(a) Thermodynamic data at 298.15 K 


Heat capacity, C p , 
cal/mole-K 

27.605 

26.719 

7.166 

Enthalpy, AjH, 
kcal/mole 

2.820 

-8.793 

9.318 

Entropy, S°, 
cal/mole-K 

83.302 

80.289 

43.881 


Heat capacity C and entropy S° at various temperatures, cal/mole-K 


Temperature 

T, 

K 

C P 

400 

36.188 

600 

49.870 

1000 

65.214 

1200 

69.204 

1400 

72.434 

1600 

75.006 

1800 

77.020 

2000 

78.563 


92.635 

110.05 

139.56 

151.82 

162.74 
172.58 
181.54 

189.74 


89.323 

106.04 

133.89 

145.31 

155.43 

164.54 

172.81 

180.37 


5 




































TABLE II.-— BENZENE OXIDATION SUBMECHANISM 


Number 

Reaction 

Forward rate coefficients 

Reference 



3 Aj ' 

n i 






cm , mole, s 


_J 

cal/mole 



1 

c 6 h 6 + o 2 -c 6 h 5 o + oh 

4.0x10° 

0 

34 000 

Bittker (1991) 


2 

C 6 H 6 + C 6 H 5-C 12 H 10 + H 

4.0x10' 1 



4 000 

Fahr and Stein (1988) 


3 

c 6 h 6 -c 6 h 5 + h 

5.0xl0 15 



108 000 

Hsu, Lin, and Lin (1984) a 

4 

c 6 h 6 + h^c 6 h 5 + h 2 

2.5xl0 14 



16 000 

Kiefer et al. (1985) 


5 

C 6 H6 + 0 ** C 6 H 5 0 + H 

2.78x10° 



4 910 

Nicovich, Gump, and Ravishankara (1982) 

6 

+ OH C 6 H 5 + H 2 0 

2.13x10° 



4 580 

Madronich and Felder (1985) 

7 

c 6 h 6 + c 2 h^c 6 h 5 + c 2 h 2 

1.0x10° 



0 

Estimated 


8 

C6H6 + c 2 H3-*c 6 h 5 + c 2 h 4 

1.0x10° 



0 

Estimated 


9 

c 6 h 6 + c 4 h 5 -c 6 h 5 + c 4 h 6 

1.0x10° 



10 000 

Estimated 


10 

c 6 h 6 + ch 3 ~c 6 h 5 + ch 4 

4.365x10^ 

5.0 

12 300 

Pamidimukkala et al. (1987) 

11 

c 6 h 5 + c 2 h 2 ~c 6 h 5 c 2 h + h 

3.24x10" 

0 

1 350 

Vaughn, Howard, and Longwell (1991) 

12 

c 4 h 3 + m ** c 4 h 2 + h + m 

l.OxlO 16 



60 000 

Miller et al. (1982) 


13 

c 6 h 5 o ** c 5 h 5 + CO 

2.51x10" 



43 900 

Lin and Lin (1986) 


14 

c 6 h 5 + o 2 -c 6 h 5 o + o 

2.1x10° 



7 470 

Lin and Lin (1987) 


15 

C 6 H 5 + H0 2 ^ C 6 H 5 0 + OH 

2.0x10° 



1 000 

Bittker (1991) 


16 

c 6 h 5 ^c 4 h 3 + c 2 h 2 

4.5x10° 



72 530 

Braun-Unkhoff, Frank, and Just (1988) 

17 

c 6 h 5 oh ^ c 6 h 5 o + H 

2.0xl0 16 



88 000 

Bittker (1991) 


18 

C^OH + H CgHg + OH 

2.2x10° 



7 910 

Lovell, Brezinsky, and Glassman (1989) 

19 

c 6 h 5 oh + h-c 6 h 5 o + h 2 

1.15xl0 14 



12 405 

Lovell, Brezinsky, and Glassman (1989) 

20 

C 5 H s + C 6 H 5 OH ^ C 6 H 5 0 + CjHg 

4.20x10° 

-0.82 

19 840 

Emdee, Brezinsky, and Glassman (1992) b 

21 

c 5 H6-*c 5 h 5 + h 

8.13xl0 24 

-2.98 

78 682 

Dean (1990) 


22 

C 5 H 6 + Oj ^ C 5 H 5 0 + OH 

1.0x10° 

0 

20 712 

Bittker (1991) 


23 

C 6 H 5 OH + OH -* c 6 h 5 o + h 2 o 

3.0x10° 



0 

Bittker (1991) 


24 

c 6 h 5 oh + ho 2 c 6 h 5 o + h 2 o 2 

3.0x10° 



1 500 

Bittker (1991) 


25 

c 6 h s oh + 0 c 6 h 5 o + OH 

2.81x10° 



7 352 

Emdee, Brezinsky, and Glassman (1992) 

26 

c 6 h 5 oh + c 2 h 3 ^c 6 h 5 o + c 2 h 4 

6.0x10° 



0 



27 

c 6 h 5 oh + c 4 h 5 ^ c 6 h 5 o + c 4 h 6 

6.0x10° 



0 



28 

c 6 h 5 oh + c 6 h 5 ^c 6 h 5 o + c 6 h 6 

4.91x10° 



4 400 



29 

C 5 H 6 + OH ** C 5 H 5 + H 2 0 

3.43xl0 9 

1.18 

-447 



30 

c 5 h 6 + h^c 5 h 5 + h 2 

2.19xl0 8 

1.77 

3 000 



31 

c 5 h 6 + o ~ c 5 h 5 + OH 

1.81x10° 

0 

3 080 



32 

c s h 6 + c 2 h 3 -c 5 h 5 + c 2 h 4 

6.0x10° 



0 



33 

C 5 H 6 + C 4 H 5-C s H 5 + C 4 H 6 

6.0x10° 



0 



34 

c 4 h 6 -c 4 h 5 + h 

1.2xl0 16 



109 910 

Vaughn, Howard, and Longwell (1991) 

35 

c 4 h 6 + oh ~ c 4 h 5 + h 2 o 

4.79x10° 



1 230 

Vaughn, Howard, and Longwell (1991) 

36 

c 4 h 6 + h~c 4 h 5 + h 2 

1.51xl0 14 



10 200 

Vaughn, Howard, and Longwell (1991) 

37 

c 5 h 6 + ho 2 -c 5 h 5 + h 2 o 2 

1.99x10° 



11 600 

Emdee, Brezinsky, and Glassman (1992) 

38 

c 5 h 5 o c 4 h 5 + CO 

3.0xl0 16 



15 000 

Bittker (1991) 


39 

c 5 h 5 + 0 - c 4 h 5 + CO 

l.OxlO 14 



0 

Emdee, Brezinsky, and Glassman (1992) 

40 

C 5 H 5 + OH C 5 H 4 OH + H 

1.0x10° 



0 

Bittker (1991) 


41 

C 5 H 4 OH ^ C 4 H 4 + HCO 

1.0x10° 



22 000 

Bittker (1991) 


42 

C 4 H 4 + C 6 H 5 -C 6 H 6 + C 4 H 3 

1.0x10° 



0 

Colket (1986) 


43 

C 4 H 4 ~ 2 C 2 H 2 

1.3x10° 



82 500 

Kiefer et al. (1988) 


44 

C 4 H 4 + O ^ C 4 H 3 + OH 

1.0x10° 



0 

Estimated 


45 

c 4 h 4 + oh ~ c 4 h 3 + h 2 o 

1.0x10° 



0 

Estimated 


46 

c 4 h 4 + h ~ c 4 h 3 + h 2 

1.6xl0 14 



14 500 

Colket (1986) 


47 

c 4 h 4 + c 2 h - c 4 h 3 + c 2 h 2 

4.0x10° 



0 

Frenklach et al. (1983) 

48 

c 5 h 5 + ho 2 c 5 h 5 o + OH 

2.0x10° 



0 

Bittker (1991) 



“Adjustment of A-factor to l.OxlO 16 for increased pressure (approx) from Kiefer et al. (1985) used for ignition delay time computations. 


b See discussion of this rate coefficient in the section Benzene Oxidation. 
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TABLE II.— Continued. BENZENE OXIDATION SUBMECHANISM 


Number 

Reaction 

Forward rate coefficients 


Reference 

Aj, 

cm 3 , mole, s 

n j 

E J’ 

cal/mole 


49 

2C 6 H s -C 12 H 10 

3.1x10° 

0 


0 

Bittker (1991) 


50 

C 4 H 5 ** C 2 H 3 + C 2 H 2 

2.0x10' 1 

0.7 

42 260 

(c) 


51 

C 4 H 5 + M -H. c 4 H 4 + H + M 

3.0xl0 15 

0 

32 000 

Est. from reaction (60) 

52 

c 4 h 2 + 0 c 2 ho + c 2 h 

l.OxlO 13 




0 

McLain, Jachimowski, and Wilson (1979) 

53 

C 4 H 2 + OH ^ C 3 H 2 + HCO 

3.0xl0' 3 




0 

Miller et al. (1982) 

54 

C 4 H 2 + O ** CO + c 3 h 2 

1.2xl0 12 




0 

Miller et al. (1982) 

55 

C 2 H 4 + M ~ C 2 H 2 + H 2 + M 

9.33xl0 16 



77 200 

Miller et al. (1982) 

56 

C 2 H 4 + OH -H. C 2 H 3 - 1 - H 2 0 

4.79x10° 



1 230 

Westbrook and Dryer (1984) 

57 

C 2 H 4 + O CH 3 + HCO 

3.31xl0 12 



1 130 



58 

c 2 h 4 + 0 ** ch 2 o + ch 2 

2.51xl0 13 



5 000 



59 

c 2 h 4 + oh ** ch 3 + ch 2 o 

2.0xl0 12 



960 



60 

c 2 h 3 + m ^ c 2 h 2 + h + m 

3.0xl0 15 



32 000 

Wamatz (1984) 


61 

C 2 H 3 + 0 2 ch 2 0 + HCO 

3.98xl0 12 



-250 

Slagle et al. (1984) 

62 

c 2 h 3 + h c 2 h 2 + H 2 

6.0xl0 12 




0 

Miller et al. (1982) 

63 

C 2 H 3 + OH ^ C 2 H 2 + H 2 0 

5.00xl0 12 







64 

c 2 h 3 + ch 2 ~ c 2 h 2 + ch 3 

3.00xl0 13 







65 

C 2 H 3 + C 2 H 2 C 2 H 2 

3.00x10° 







66 

C 2 H 3 + O -H. C 2 H 2 0 + H 

3.30x10° 







67 

ch 2 + ch 2 ** c 2 h 2 + h 2 

4.00x10° 







68 

CH 2 + CH 2 ~ C 2 H 3 + H 

5.00x10° 





Westbrook and Dryer (1984) 

69 

CH 2 + OH ** CH + H 2 0 

2.51x10" 

0.67 

25 700 

Miller et al. (1982) 

70 

CH 2 + O ~ CH + OH 

2.00x10" 

.68 

25 000 



71 

CH + 0 2 - C0 2 + 2H 

1.59x10° 

0 

1 000 



72 

C 2 H 2 + M -h. c 2 h + h + m 

4.17X10 16 



107 000 



73 

c 2 h 2 + c 2 h 2 c 4 h 3 + H 

2.00x10° 



45 900 



74 

C 2 H 2 + C 2 H ** C 4 H 2 + H 

3.00x10° 




0 



75 

c 2 h 2 + O ch 2 + CO 

1.60xl0 14 



9 890 

Frank, Bhaskaran, and Just (1986) 

76 

c 2 h 2 + 0 ~ c 2 ho + H 

4.00xl0 14 



10 660 

Frank, Bhaskaran, and Just (1986) 

77 

c 2 h 2 + oh c 2 h + h 2 o 

6.31x10° 



7 000 

Westbrook and Dryer (1984) 

78 

C 2 H 2 + OH ^ C 2 H 2 0 + H 

3.20x10" 



200 

Westbrook and Dryer (1984) 

79 

C 2 H 2 + CH 2 ^ C 3 H 3 + H 

1.20x10° 



6 600 

Boehland, Temps, and Wagner (1986) 

80 

C 3 H 4 + M ~ C 3 H 3 + H + M 

2.00x10° 



65 000 

Pamidimukkala et al. (1987) 

81 

C 3 H 3 + 0 2 ^ C 2 H 2 0 + HCO 

3.00x10'° 



2 870 

Slagle and Gutman (1986) 

82 

C 2 H 2 0 + OH ^ CH 2 0 + HCO 

2.80x10° 




0 

Miller et al. (1982) 

83 

C 2 H 2 0 + OH ** C 2 HO + H 2 0 

7.50x10° 



3 000 

Miller et al. (1982) 

84 

c 2 h 2 o + H ~ CH 3 + CO 

1.13x10° 



3 428 

Westbrook and Dryer (1984) 

85 

c 2 h 2 o + h c 2 ho + h 2 

7.50x10° 



8 000 

Miller et al. (1982) 

86 

C 2 H 2 0 + O C 2 HO + OH 

5.00x10° 



8 000 



87 

c 2 h 2 o + O - ch 2 o + CO 

2.00x10° 



0 



88 

c 2 h 2 o + m ^ ch 2 + CO + m 

2.00X10 16 



60 000 



89 

C 2 HO + 0 2 2 CO + OH 

1.46x10° 



2 500 



90 

C 2 HO + O -H. 2 CO + H 

1.20x10° 



0 

Westbrook and Dryer (1984) 

91 

C 2 HO + OH ^ 2 HCO 

1.00x10° 



0 

Miller et al. (1982) 

92 

C 2 HO + H 1 * CH 2 + CO 

5.00x10° 



0 



93 

c 2 ho + ch 2 ^ c 2 h 3 + CO 

3.00x10° 



0 



94 

c 2 ho + ch 2 ^ ch 2 o + c 2 h 

1.00x10° 



2 000 



95 

2 C 2 HO ^ C 2 H 2 + 2 CO 

1.00x10° 



0 



96 | 

c 2 h + oh ** c 2 ho + H 

2.00x10° 



0 




c Modified from Emdee, Brezinsky, and Giassman (1992) by reducing ,4 -factor. 
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TABLE II.— Continued. BENZENE OXIDATION SUBMECHANISM 


Number 

Reaction 

Forward rate coefficients 

Reference 




n 

/ 


H 





cm 3 , mole, s 


cal/mole 



97 

C 2 H + 0 2 '** C 2 HO + O 

5.00xl0 13 

0 

1 500 

Miller et al. (1982) 

98 

c 2 h + o **co + ch 

5.00xl0 13 




0 

Miller et al. (1982) 

99 

ch 4 + m ^ ch 3 + H + m 

2.00xl0 !7 



88 000 

Westbrook and Dryer (1984) 

100 

ch 4 + o 2 ch 3 + ho 2 

7.94xl0 13 



56 000 

Westbrook and Dryer (1984) 

101 

CH 4 + H ** CH 3 + H 2 

1.26xl0 14 



11 900 

Brabbs and Brokaw (1974) 

102 

ch 4 oh - ch 3 + h 2 o 

2.50xl0 13 



5 010 

Brabbs and Brokaw (1974) 

103 

ch 4 + o ~ ch 3 + oh 

1.90x10 14 



11 720 

Brabbs and Brokaw (1974) 

104 

CH 3 + 0 2 -h. ch 3 o + o 

4.79xl0 13 



29 000 

Westbrook and Dryer (1984) 

105 

ch 3 + oh ** ch 3 o + H 

6.30xl0 12 




0 

Westley (1980) 


106 

CHjO + M ^ CH 2 0 + H + M 

5.00xl0 13 



21 000 

Westbrook and Dryer (1984) 

107 

ch 3 o + o 2 ** ch 2 o + ho 2 

l.OOxlO 12 



6000 



108 

CHjO + H ** CH 2 0 + H 2 

2.00xl0 13 




0 



109 

ch 3 + ch 3 C 2 H 4 + h 2 

l.OOxlO 16 



32 000 



110 

CH 3 + O ^ CH 2 0 + H 

1.29xl0 14 



2 000 



111 

CH 3 + CH 2 0 ** CH 4 + HCO 

l.OOxlO 10 

0.5 

6 000 



112 

CH 3 + HCO - CH 4 + CO 

3.00X10 11 

0.5 


0 



113 

CH 3 + H0 2 CHjO + OH 

2.00xl0 13 

0 


0 



114 

CH 2 0 + M H + HCO + M 

S.OOxlO 16 



76 500 

Wamatz (1984) 

115 

CH 2 0 + OH HCO + H 2 0 

3.00xl0 13 



1 200 



116 

CH 2 0 + H - HCO + H 2 

2.5Qxl0 13 



3 990 



117 

CH 2 0 + O ** HCO + OH 

3.50xl0 13 



3 510 



118 

HCO + H0 2 ** CH 2 0 + 0 2 

l.OOxlO 14 



3 000 

Westley (1980) 


119 

HCO + M-*H + CO + M 

2.94X10 14 



15 569 

Cherian et al. (1981 ) d 

120 

HCO + 0 2 - CO + H0 2 

3.31x10° 



7 000 

Westbrook and Dryer (1984) 

121 

hco + oh**co + h 2 o 

l.OOxlO 14 




0 



122 

HCO + H CO + H 2 

2.00xl0 14 







123 

HCO + O ** CO + OH 

l.OOxlO 14 







124 

CH + 0 2 v* HCO + O 

1.00x10° 







125 

CO + O + M ** C0 2 + M 

5.90x10° 



4 100 



126 

CO 4- O2 ** CO2 4 O 

2.50x10° 



47 690 



127 

CO + OH ^ C0 2 + H 

4.17x10° 



1 000 

Brabbs and Brokaw (1974) 

128 

CO + H0 2 ** C0 2 + OH 

5.75x10° 



22 930 

Westbrook and Dryer (1984) 

129 

O + H 2 0 ** OH + OH 

6.80x10° 



18 365 

Brabbs and Musiak (1988) 

130 

H + 0 2 v. OH + O 

1.89xl0 14 



16 400 



131 

O + H 2 OH + H 

4.20xl0 14 



13 750 



132 

H + H0 2 -h> H 2 + 0 2 

7.28x10° 



2 126 



133 

O + H0 2 ** OH + 0 2 

5.00x10° 



1 000 



134 

H0 2 + OH H 2 0 + 0 2 

8.00x10° 




0 



135 

H + H0 2 OH + OH 

1.34xl0 14 



1 070 




d Computed from reverse rate coefficient and equilibrium constant. 


Collisionat Efficiencies 


Reaction (140): H 2 = 2.3; 0 2 = 0.78; H 2 0 = 6.0; H 2 0 2 = 6.6 

Reaction (142): H 2 = 3.0; 0 2 = 1.3; H 2 0 = 21.3; N 2 = 1.3; C0 2 = 7.0; C 6 H 6 = 20.0; CH 4 = 5.0 

Reaction (143): H 2 = 4.0; 0 2 = 1.5; H z O = 20.0; N 2 = 1.5; C0 2 = 4.0; C 6 H 6 = 20.0 

Reaction (145): H 2 = 4.1; 0 2 = 2.0; H 2 0 = 15.0; N 2 - 2.0 




TABLE II.— Concluded. BENZENE OXIDATION SUBMECHANISM 


Number 

Reaction 

Forward rate coefficients 

Reference 

4 h 

cm 3 , mole, s 

n. 

j 

E r 

cal/mole 

136 

h 2 + ho 2 ^h 2 o 2 + h 

7.91xl0 13 

0 

25 000 

Brabbs and Musiak (1988) 

137 

OH + H 2 0 2 H 2 0 + H0 2 

6-lOxlO 12 



1 430 



138 

ho 2 + ho 2 ** h 2 o 2 + o 2 

1.80xl0 12 



0 



139 

H + H 2 0 2 -h. OH + H 2 0 

7.80x10" 



0 



140 

H 2 0 2 + M ^ OH + OH + M 

1.44xl0 17 



45 510 



141 

H 2 + OH H 2 0 + H 

4.74xl0 13 



6 098 



142 

H + 0 2 + M H0 2 + M 

1.46xl0 15 



-1 000 



143 

h 2 o + M^H + OH + M 

1.30xl0 15 



105 140 



144 

O + H + M^OH + M 

7.10xl0 18 

-1.0 

0 



145 

h 2 + m-*h + h + m 

2.20x1 0 14 

0 

96 000 



146 

O 2 + M O + O + M 

1.80xl0 18 

-1.0 

118 020 




TABLE III.— INITIAL CONDITIONS FOR OXIDATION OF BENZENE AND TOLUENE 
IN FLOW REACTOR AT 101 kPa (1 atm) 


Case 

number 

Fuel 

Mixture concentrations, 
ppmv 

Temperature, 

T, 

K 

Equivalence 

ratio, 21 

<P 

Fuel 

Oxygen 

Nitrogen 

B-l 

Benezene b 

1571 

15 900 

982 529 

1098 

0.74 

B-2 

Benzene 0 

1495 

12 301 

986 204 

1096 

0.91 

B-3 

Benzene b 

1591 

11 700 

986 709 

1096 

1.0 

B-4 

Benzene b 

1581 

8 700 

989 719 

1096 

1.36 

T-l 

Toluene 0 

1540 

20 015 

978 445 

1188 

0.69 

T-2 

Toluene 0 

1616 

10 907 

987 477 

1190 

1.33 


a (Fuel -oxygen mole ratio)/(stoichiometric fuel-oxygen mole ratio). 
Novell, Brezinsky, and Glassman (1988). 
c Emdee, Brezinsky, and Glassman (1992). 





TABLE IV.— COMPARISON OF COMPUTED AND EXPERIMENTAL IGNITION 
DELAY TIMES FOR BENZENE OXIDATION BEHIND A REFLECTED SHOCK 


Mixture number 
and description 

Initial 

temperature, 

T, 

K 

Ignition delay time, t, 
ps 

Percentage 

difference 

Experimental 

Computed 

1 

1209 

878 

870 

-0.9 

Equivalence ratio, <p, 0.5 

1227 

743 

721 

-3.0 

Mole % argon, 78.333 

1254 

435 

524 

20.5 

Pressure, p, kPa (atm) 

1276 

330 

416 

26.1 

= 190.9 to 231.4 

1291 

272 

360 

32.4 

(1.890 to 2.291) 

1314 

202 

282 

39.6 


1345 

159 

210 

32.1 


— 

— 

... 

a 26.1 

2 

1345 

755 

750 

-0.7 

Equivalence ratio, <p, 1.0 

1374 

604 

570 

-5.6 

Mole % argon, 95.616 

1402 

415 

450 

8.4 

Pressure, p, kPa (atm) 

1412 

412 

413 

0.2 

= 573.1 to 720.5 

1428 

367 

358 

-2.4 

(5.674 to 7.134) 

1482 

213 

230 

8.0 


1525 

122 

161 

32.0 


1528 

122 

158 

29.5 


— 

— 

... 

a l 6. 1 

3 

1283 

750 

530 

-29.3 

Equivalence ratio, <p, 1 .0 

1290 

613 

496 

-19.1 

Mole % argon, 85.635 

1294 

607 

466 

-23.2 

Pressure, p, kPa (atm) 

1328 

490 

330 

-32.7 

= 205.5 to 250.5 

1355 

303 

260 

-14.2 

(2.035 to 2.480) 

1369 

287 

232 

-19.4 


1379 

291 

211 

-27.5 


1405 

198 

170 

-14.1 


1408 

189 

164 

-13.2 


1417 

178 

153 

-14.0 


1435 

151 

127 

-15.9 


— 

— 

-- 

*21.3 

4 

1363 

1520 

690 

-54.6 

Equivalence ratio, <p, 2.0 

1415 

890 

458 

-48.5 

Mole % argon, 93.553 

1457 

599 

312 

-47.9 

Pressure, p, kPa (atm) 

1540 

274 

168 

-38.7 

= 202.6 to 265.3 

1554 

243 

150 

-38.3 

(2.006 to 2.627) 

1570 

211 

132 

-37.4 


1582 

157 

123 

-21.7 


1600 

154 

107 

-30.5 


— 

™ 

... 

a 40.9 


‘Standard deviation, o. 




TABLE V.— TEMPERATURE DEPENDENCE OF IGNITION 
DELAY TIMES FOR BENZENE OXIDATION BEHIND A 
REFLECTED SHOCK 


Mixture number and description 

Activation energy factor, 
AE, cal 

Experimental 

Computed 

1 

Equivalence ratio, <p, 0.5 
Pressure, p, kPa (atm) = 202 (2) 

42 718 

33 949 

2 

Equivalence ratio, tp, 1.0 (dilute) 
Pressure, p, kPa (atm) = 606 (6) 

41 472 

34 828 

3 

Equivalence ratio, <p, 1 (strong) 
Pressure p, kPa (atm) = 202 (2) 

37 251 

33 649 

4 

Equivalence ratio, <p, 2.0 
Pressure, p, kPa (atm) = 202 (2) 

42 404 

34 264 




TABLE VI.— TOLUENE OXIDATION SUBMECHANISM 


Number 

Reaction 

Coefficients in the equation 

Reference 



k = AjT 

" y exp (Ej/RT) 




A j’ 

n 


E r 




cm 3 , mole, s 



cal/mole 


■1 

c 6 h 5 ch 3 + o 2 c 6 h 5 ch 2 + ho 2 

2.50xl0 14 


41 400 

This work 

IB 

h + c 6 h 5 ch 2 ^c 6 h 5 ch 3 

1.80xl0 14 



0 

Emdee, Brezinsky, and Glassman (1992) 

IB 

c 6 h 5 ch 3 -c 6 h 5 + ch 3 

1.40xl0 16 



99 800 



IB 

c 6 h 5 + c 6 h 5 ch 3 - c 6 h 6 + c 6 h 5 ch 2 

2.10xl0 12 



4 400 



M5M 

ch 3 + c 6 h 5 ch 3 ~c 6 h 5 ch 2 + ch 4 

3.16xl0 n 



9 500 



152 

c 6 h 5 ch 3 +h-c 6 h 6 + ch 3 

1.50xl0 13 



5 148 

This work 

153 

c 6 h 5 ch 3 + h ^ c 6 h 5 ch 2 + h 2 

l.OOxlO 14 



8 235 

This work 

154 

c 6 h 5 ch 3 + OH ** c 6 h 5 ch 2 + h 2 o 

1.26xl0 13 



2 583 

Emdee, Brezinsky, and Glassman (1992) 

155 

c 6 h 5 ch 3 + o ch 3 c 6 h 4 o + H 

1.63xl0 13 



3 418 

Emdee, Brezinsky, and Glassman (1992) 

156 

ch 3 c 6 h 4 o + h ch 3 c 6 h 4 oh 

2.50xl0 14 



0 

Emdee, Brezinsky, and Glassman (1992) 

157 

ch 3 c 6 h 4 o ~ c 6 h 6 + H + CO 

6.50xl0 10 



43 900 

This work 

158 

ch 3 c 6 h 4 oh + oh ch 3 c 6 h 4 o + h 2 o 

6.00X10 12 



0 

Emdee, Brezinsky, and Glassman (1992) 

159 

ch 3 c 6 h 4 oh + H ^ ch 3 c 6 h 4 o + h 2 

1.15xl0 14 



12 400 



160 

ch 3 c 6 h 4 oh + h - c 6 h 5 ch 3 + OH 

2.21xl0 13 



7 910 



161 

ch 3 c 6 h 4 oh + H ** c 6 h 5 oh + ch 3 

1.20xl0 13 



5 148 

I 


162 

CH,C fi H 4 OH + C fi H s CH, ** CH,QH 4 0 + C fi H s CH, 

1.05xl0 n 



9 500 



163 

c 6 h 5 ch 2 + o ^ c 6 h 5 CHO + H 

3.75xl0 14 



0 

This work 

164 

c 6 h 5 ch 2 + o~c 6 h 5 + ch 2 o 

8.00x10 !3 



0 

Emdee, Brezinsky, and Glassman (1992) 

165 

c 6 h 5 ch 2 + oh ~ c 6 h 5 ch 2 oh 

6.00xl0 13 



0 



166 

C 6 H 5 CH 2 OH + 0 2 -Mi c 6 h 5 cho + h + ho 2 

2.00xl0 14 



41 400 



167 

C 6 H 5 CH 2 OH + H ** C 6 H6 + ch 2 oh 

1.20xl0 13 



5 148 



168 

c 6 h 5 ch 2 oh + OH ** C 6 H 5 CHO + H + h 2 o 

8.43xl0 12 



2 583 



169 

c 6 h 5 ch 2 oh + h - c 6 h 5 cho + h + h 2 

8.00x10° 



8 235 



170 

C fi H s CH ? OH + C fi H s CH, C fi H,CHO + C fi H s CH, + H 

2-llxlO 11 



9 500 



171 

c 6 h 5 ch 2 oh + c 6 h 5 c 6 h 5 cho + c 6 h 6 + H 

1.40xl0 12 



4 400 



172 

CH 3 + OH ** CH 2 OH + H 

1.09xl0 n 

0.40 

-708 



173 

ch 2 oh + o 2 - ch 2 o + ho 2 

2.41xl0 14 

0 

5 000 



174 

CH 2 OH + M CH 2 0 + H + M 

1.67xl0 24 

-2.5 

34 190 



175 

CH 2 0 + H0 2 v. HCO + H 2 0 2 

1.99xl0 12 

0 

11 600 



176 

c 6 h s ch 2 + ho 2 c 6 h 5 cho + H + OH 

2.50xl0 14 

0 

0 



177 

2 C 6 H 5 CH 2 bibenzyl 

2.51xl0 n 

0.4 

0 



178 

C 6 H 5 C 2 H 5 ** c 6 h 5 ch 2 + ch 3 

2-OOxlO 15 

0 

72 700 



179 

C fi H s C,H s + OH ** C fi H s C,H, + H + H,0 

8.43xl0 12 



2 583 



180 

c 6 h 5 c 2 h 5 + H- c 6 h 5 c 2 h 3 + H + h 2 

8.00x10° 



8 235 



181 

c 6 h 5 c 2 h 5 + o 2 c 6 h 5 c 2 h 3 + h + ho 2 

2.00x10° 



41 400 



182 

c 6 h 5 cho + o 2 ~ c 6 h 5 co + ho 2 

1.02x10° 



38 950 



183 

c 6 h 5 cho + oh c 6 h 5 co + h 2 o 

1.71xl0 9 

1.18 

-447 



184 

c 6 h 5 cho + h ** c 6 h 5 co + h 2 

5.00x10° 

0 

4 928 



185 

C 6 H 5 CHO + H -h. c 6 h 6 + HCO 

1.20x10° 

0 

5 148 



186 

c 6 h 5 cho + o ~ c 6 h 5 co + OH 

9.04xl0 12 

0 

3 080 



187 

C 6 H 5 CHO + c 6 h s ch 2 ** c 6 h 5 co + c 6 h 5 ch 3 

2.77xl0 3 

2.81 

5 773 



188 

c 6 h 5 cho + ch 3 c 6 h 5 co + ch 4 

2.77xl0 3 

2.81 

5 773 



189 

c 6 h 5 cho + c 6 h 5 ** c 6 h 5 co + c 6 h 6 

7-OlxlO 1 1 

0 

4 400 



190 

c 6 h 5 co-c 6 h 5 + co 

4.00xl0 14 

0 

29 400 



191 

c 6 h 5 oh + c 6 h 5 ch 2 c 6 h 5 o + c 6 h 5 ch 3 

1.05x10° 

0 

9 500 
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TABLE VII.— COMPARISON OF COMPUTED AND EXPERIMENTAL IGNITION DELAY 


TIMES FOR TOLUENE OXIDATION BEHIND A REFLECTED SHOCK 


Mixture number and 
description 

Initial 

temperature, 

T, 

K 

Ignition delay time, x 

]iS 

Percentage 

difference 

Experimental 

Computed 

1 

1334 

847 

975 

15.1 

Equivalence ratio, <p, 0.33 1 

1353 

657 

760 

15.7 

Mole % argon, 85.989 

1358 

556 

702 

26.3 

Pressure, p, kPa (atm), 

1385 

388 

493 

27.1 

197.0 to 236.3 (1.95 to 2.34) 

1406 

287 

383 

33.4 


1419 

221 

320 

44.8 


1432 

186 

276 

48.4 


1437 

183 

255 

39.3 


1443 

124 

228 

83.9 


— - 

— 

... 

a 42.1 

2 

1426 

1311 

1140 

-13.0 

Equivalence ratio, <p, 1 .0 

1436 

1104 

1019 

-7.7 

Mole % argon, 95.027 

1441 

990 

948 

-4.2 

Pressure, p, kPa (atm), 

1496 

452 

475 

5.1 

198.0 to 241.4 (1.96 to 2.39) 

1505 

452 

426 

-5.8 


1505 

370 

430 

16.2 


1520 

335 

364 

8.7 


1561 

207 

218 

5.3 


1611 

101 

120 

18.8 


— 

— 

— 

a 10.7 

3 

1382 

741 

870 

17.4 

Equivalence ratio, (p 1.0 

1417 

563 

580 

3.0 

Mole % argon, 95.027 

1433 

490 

488 

-0.4 

Pressure, p, kPa (atm), 

1476 

330 

300 

-9.1 

566.6 to 674.7 (5.61 to 6.68) 

1516 

237 

190 

-19.8 


1520 

208 

183 

-12.0 


1540 

175 

146 

-16.6 


.... 

— 

— 

“13.1 

4 

1353 

1030 

582 

-43.5 

Equivalence ratio, cp 1.0 

1362 

975 

520 

-46.7 

Mole % argon, 85.053 

1390 

596 

370 

-37.9 

Pressure, p, kPa (atm), 

1423 

356 

259 

-27.2 

236.3 to 293.9 (2.34 to 2.91) 

1442 

328 

208 

-36.6 


1483 

185 

128 

-30.8 


1535 

100 

73 

-27.0 


.... 

— 

— 

“36.4 


‘Standard deviation, a. 




TABLE Vin.— TEMPERATURE DEPENDENCE OF 
IGNITION DELAY TIMES FOR TOLUENE 
OXIDATION BEHIND A REFLECTED SHOCK 


Mixture number 
and description 

Activation energy factor, A E, cal 

Experimental 

data 

Computed 

results 

1 

Equivalence ratio, q>, 0.331 
Pressure, p, kPa (atm), = 202 (2) 

61 853 

50 043 

2 

Equivalence ratio, q> 1.0 (dilute) 
Pressure, p, kPa (atm), = 202 (2) 

61 774 

55 258 

3 

Equivalence ratio, <p 1.0 (dilute) 
Pressure, p, kPa (atm), = 606 (6) 

38 786 

47 872 

4 

Equivalence ratio, <p 2.0 (strong) 
Pressure, p, kPa (atm), = 202 (2) 

53 256 

46 864 




Species concentration, ppmv Species concentration, ppmv Species concentration, ppmv 



Experiment 

Computed 

New mechanism 

Original mechanism 



Experiment 

Computed 

New mechanism 

Original mechanism 



Experiment 

Computed 
New mechanism 
Original mechanism 


Time, ms 

Figure 1 . — Concentration time profiles of benzene and carbon monoxide for benzene oxidation 
in flow reactor, (a) cp = 0.74; T 0 = 1 098 K. (b) <p = 1 .0; T 0 = 1 096 K. (c) cp = 1 .36; T 0 = 1 096 K. 



Species concentration, ppmv Species concentration, ppmv Species concentration, ppmv 


240 r— 



Experiment 

Computed 
New mechanism 
Original mechanism 



Experiment 

Computed 
New mechanism 
Original mechanism 



Experiment 

Computed 
New mechanism 


Figure 2. — Concentration time profiles of cyclopentadiene and phenol for benzene oxidation 
in flow reactor, (a) <p = 0.74; T 0 = 1 098 K. (b) <p = 1.0; T 0 = 1 096 K. (c) <p - 1 .36; T 0 = 1 096 K. 
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Figure 3. — Concentration time profiles of benzene and phenol for benzene 
oxidation in flow reactor; <p = 0.91 ; T 0 = 1 096 K. 



Figure 4. — Concentration time profiles of carbon monoxide and C 2 hydrocarbons for 
benzene oxidation in flow reactor; (p = 0.91 ; T 0 = 1 096 K. 



Ignition delay time, t, p,s 


Experiment 



Experiment 


Computed 



Time, ms 


Figure 5. — Concentration time profiles of cyclopentadiene and C 4 hydrocarbons for 
benzene oxidation in flow reactor; (p = 0.91 ; T 0 = 1 096 K. 



1350 1280 1219 

Temperature, T, K 

Figure 6.— Ignition delay times for benzene oxidation behind reflected shock; mixture 1 ; 
<p = 0.5; mole % argon = 78.333. 






























































































































Species concentration, ppmv Species Species concentration, ppmv 

concentration, ppmv 




Experiment 

Computed 
This work 
Emdee 


Experiment 

Computed 

This work 

Emdee 



Experiment 


Computed 

This work 

— Emdee 


Figure 1 2. — Concentration time profiles of toluene and benzene for 
toluene oxidation in flow reactor, (a) and (b) <p = 0.69; T 0 = 1 1 88 K. 
(c) <p = 1.33; T 0 = 1190 K. 



160 


Experiment 



Figure 1 3. — Concentration time profiles of benzaldehyde and methane for 
toluene oxidation in flow reactor; <p = 0.69; T 0 = 1 188 K. 



Figure 14. — Concentration time profile of benzyl alcohol for toluene oxi- 
dation in flow reactor; = 0.69; T 0 = 1 188 K. 
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Figure 15. — Concentration time profiles of benzaldehyde and benzyl alcohol for 
toluene oxidation in flow reactor; <p = 1 .33; T 0 = 1 190 K. 



Figure 1 6. — Concentration time profiles of carbon monoxide and acetylene 
for toluene oxidation in flow reactor; <p = 0.69; T 0 = 1188 K. 






Cyclopentadiene concentration, ppmv Cresols concentration, ppmv 




Figure 18. — Concentration time profiles of cresols and cyclopentadiene for 
toluene oxidation in flow reactor; <p = 0.69; T 0 = 1 188 K. 











Ignition delay time, ps 



10 000/T, K" 1 


1610 1515 1430 1350 

Temperature, T, K 

Figure 21 . — Ignition delay times for toluene oxidation behind reflected shock; 
mixture 3; <p = 1.0; mole % argon = 95.207; pressure = 606 kPa (= 6 atm). 



1610 1515 1430 1350 

Temperature, T, K 

Figure 22. — Ignition delay times for toluene oxidation behind reflected 
shock; mixture 4; <p = 1 .0; mole % argon = 85.053; pressure = 202 kPa 
(= 2 atm). 
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I. (c) Cresol. (d) Benzyl alcohol. 
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